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1.1 Introduction 

Most recent nucleosynthesis parameter studies ||, 11] place questions on the ability of 
high entropy neutrino wind scenarios in type II supernovae to produce r-process nuclei for 
A < 110 in correct amounts. In addition, it remains an open question whether the entropies 
required for the nuclei with A > 110 can actually be attained in type II supernova events. 
Thus, an alternative or supplementary r-process environment is needed, leading possibly to 
two different production sites for r-process nuclei: a high entropy, high Y e (neutrino wind in 
type II supernovae) and a low entropy, low Y e (decompression of neutron star (ns) material) 
scenario. 

Further indications for a production site possibly different from SN II arise from observations 
of low metallicity stars ||. It seems that the production of r-process nuclei is delayed in 
comparison with the major SN II yields, a fact that would be consistent with the merging 
scenario of two neutron stars. 

The tidal disruption of a ns by a black hole and possible consequences for nucleosynthesis 
has first been studied by Lattimer and Schramm the merging of a neutron star binary 

has been discussed by Symbalisty and Schramm fig] . The related decompression of the 
neutron star material has been investigated by Lattimer et al. Eichler et al. 0], who also 



discussed various other aspects of such a merging scenario, and by Meyer |1C]. In the context 



of numerical simulations the merging event nucleosynthesis has been discussed by Davies et 



al. ffl and by Ruffert et al. M]. 



1.2 The Calculations 

To investigate the possible relevance of neutron star mergers for the r-process nucleosynthe- 
sis we perform 3D Newtonian SPH calculations of the hydrodynamics of equal (1.6 M of 
baryonic) mass neutron star binary coalescences. Starting with an initial separation of 45 
km we follow the evolution of matter for 12.9 ms. We use the physical equation of state of 
Lattimer and Swesty |8| to model the microphysics of the hot neutron star matter. To test 
the sensitivity of our results to the chosen approaches and approximations we perform in total 
10 different runs where we test each time the sensitivity to one property of our model ||14|| . 
We vary the resolution (~ 21000 and ~ 50000 particles), the equation of state (polytrope), 
the artificial viscosity scheme [^], the stellar masses (1.4 M Q of baryonic matter), we include 
neutrinos (free-streaming limit), switch off the gravitational backreaction force, and vary the 
initial stellar spins. In addition we test the influence of the initial configuration, i.e. spherical 
stars versus corotating equilibrium configurations. 
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1.3 Results 



We find that, dependent on the initial spins and strongly dependent on the EOS, between 
4-10 -3 and 4- 10 -2 M become unbound. Assuming a core collapse supernova rate of 2.2 -10 -2 
(year galaxy) -1 [13|, one needs 1CT 6 to 1(T 4 M of ejected r-process material per supernova 
event to explain the observed abundances if type II supernovae are assumed to be the only 
source. The rate of neutron star mergers has recently been estimated to be 8 • 10 -6 (year 
galaxy) -1 (see [ IT | ) . Taking these numbers, one would hence need ~ 3- 10 -3 M to ~ 0.3 M 
for an explanation of the observed r-process material exclusively by neutron star mergers. 
Thus our results for the ejected mass from 4 • 10 -3 to 3 — 4 • 10 -2 M look very promising 
(see Figure II]) . 
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O SN: Ratnatunga, van den Bergh 1989, Cowan et al. 1991 

X NSM: van den Heuvel, Lorimer 1996 
* NSM: Narayan et al. 1991 
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Figure 1: The shaded region shows the amount of ejecta found in our calculations. The circle 
shows the amount of ejecta needed per event if SN II are assumed to be the only sources of 
the r-process. The asterisk gives the needed ejecta per merging event for the rate of Narayan 
et al. (1991), the cross for the estimate of van den Heuvel and Lorimer (1996). The event 
rate is given in year -1 galaxy -1 , the ejected mass in solar units. 

As a first step we use the mean properties of all ejected particles (initial corotation) for 
an r-process calculation. We adopt the following approach: in the very first expansion phase 
(where p > Pdrip ~ 4 • 10 11 g cm -3 ) we use the abundances of neutrons, protons, alphas and 
a representative nucleus provided by the LS-EOS. When the density drops below pdrip we 
switch over to a treatment of individual nuclei with a full r-process network following all reac- 
tions like neutron capture, photo-disintegrations, /3-decays etc. as discussed in Freiburghaus 
et al. Since the representative nucleus at pdrip was too neutron rich ((Z,A) = (26, 155)), 
we took the most neutron rich nucleus in the network ((Z,A) = (26,73)) and assumed the 
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remaining neutrons to be free. Following the trajectory given by the hydro calculation (ex- 
trapolation for t > 12.9 ms) we obtained the abundance pattern that is shown in Figure || 
together with the observed r-process abundances. 
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Figure 2: Comparison of the r-process calculations for a corotating system (initial corotation; 
line) with the observed abundances (crosses). 

The observed features of the abundance pattern in the range 125 < A < 200 are well 
reproduced. Especially the peak around A = 195 is easily reproduced without any tuning of 
the initial entropy. 

This approach has two shortcomings: (i) as long as the LS-EOS is used, only one (repre- 
sentative) nucleus is used instead of an ensemble of nuclei and (ii) weak interactions such as 
/3-decays or e~ — , e + — captures on protons and neutrons are disregarded in this early phase. 
For the case of initial corotation this approximation might not be crucial since the ejecta 
essentially stay cold (until perhaps, at later times, heating by /3-decays sets in). For different 
initial spins, however, weak interactions might change the Y e of the composition in this early 
phase. 

Clearly, in future investigations these aspects have to be treated in more detail. 
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